Ab initio study of ferroelectric domain walls in PbTiOa 
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We have investigated the atomistic structure of the 180° and 90° domain boundaries in the fer- 
roelectric perovskite compound PbTiOs using a first-principles ultrasoft-pseudopotential approach. 
For each case we have computed the position, thickness and creation energy of the domain walls, and 
an estimate of the barrier height for their motion has been obtained. We find both kinds of domain 
walls to be very narrow with a similar width on the order of one to two lattice constants. The energy 
of the 90° domain wall is calculated to be 35mJ/m^, about a factor of four lower than the energy 
of its 180° counterpart, and only a miniscule barrier for its motion is found. As a surprising feature 
we detected a small offset of 0.15-0.2 eV in the electrostatic potential across the 90° domain wall. 

PACS numbers: 77.80.Dj, 77.84.Dy, 61.72.Mm, 68.35.Ct 



I. INTRODUCTION 

After cooling below the transition temperature, ferro- 
electric crystals usually adopt a very complex microstruc- 
ture consisting of regions with different orientations of 
the spontaneous polarization. These domains are formed 
to reduce uncompensated depolarization fields at the sur- 
face, releasing elastic strain and so lowering the free en- 
ergy of the crystal. The domain structure and the prop- 
erties of the domain boundaries play an important role 
in the performance of many ferroelectric devices. Me- 
chanical and electrical characteristics, such as the per- 
mittivity, coercive field, and piezoelectric constants, are 
often significantly infiuenced. In particular, the thick- 
ness and the interfacial energy of the domain walls are 
important parameters in understanding the switching ki- 
netics and fatigue mechanism in ferroelectric materials. 
The width affects the wall mobility, and the energy de- 
termines how easily new domain walls may be introduced 
during polarization reversal processes. Thus, for a thor- 
ough understanding of the physical processes associated 
with the switching and fatigue behavior of a ferroelec- 
tric material, an accurate microscopic description of the 
underlying domain walls and their dynamics is required. 

A phenomenological level of description has been used 
in most previous theoretical investigations of ferroelectric 
domain walls. Several authors have modeled the struc- 
ture of domain walls-.bY using a Landau-Ginzburg type 
of continuum theory.EJ El The domain wall is described by 
a soliton solution, giving the polarization profile and the 
distribution of the elastic strain across the domain wall, 
from which the domain wall width and energy can ha 
deduced. For example, in a very early work Zhirnovlil 
obtained estimates for the thickness of domain walls in 
BaTiOa. He found the 180° domain wall to be atomi- 
cally sharp with a width of only 5-20 A, but he predicted 
the 90° domain wall to be much broader with a thick- 
ness of 50-100 A. This different picture of 180° and 90° 
domain walls has been quite commonly accepted since 



then.Q However, some of the parameters needed in the 
expansion of the free energy functional can only roughly 
be estimated from the experimental data available, the 
empirical approaches based on Landau-Ginzburg theory 
are rather limited and tend to be more qualitative in na- 
ture. Also, for atomically sharp domain walls the appli- 
cability of continuum theoretical models is questionable. 
Thus, for a deeper theoretical understanding of domain 
walls, more accurate atomistic calculations are need. 

A first work in this direction was undertaken by 
Padilla, Zhong, and VanderbiltQ who investigated 180° 
domain walls in tetragonal BaTiOs using an effective 
Hamiltonian derived from first- principles calculations. 
With this effective Hamiltonian, Monte-Carlo simula- 
tions at finite temperatures were performed to calculate 
the energy, free energy and thickness of the wall. A very 
narrow domain wall with energy 16mJ/m^ was found. 

With a first-principles pseudopotential method, based 
on density-functional theory, Poykko and ChadB stud- 
ied the 180° domain wall in PbTiOs. They confirmed a 
picture in which 180° walls are very narrow, but found a 
much higher domain- wall energy than for BaTiOs. 

In this paper, we undertake a first-principles study of 
the 180° and 90° domain walls in PbTiOa. For both 
domain wall types we have calculated the fully relaxed 
atomic structure, the domain-wall energy, the polariza- 
tion profile across the domain wall, and an estimate of 
the barrier height for their motion. As for the 180° do- 
main wall, our study is similar in spirit to the work of 
Poykko and Chadi, but the conclusions differ in some 
important respects. However, the main focus of this pa- 
per will be on the 90° domain wall, where we are not 
aware of any previous first-principles study for a ferro- 
electric perovskite. The 90° domain wall is of particular 
interest because some detailed experimental results on 
its atomistic structure have recently become available. 
(The experimental investigation of the 180° domain wall 
is much more difficult.) Some of the experiments on the 
90° domain wall have created considerable doubt regard- 
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ing the commonly accepted picture that the 90° domain 
wall is much broader than its 180° counterpart. With the 
present simultaneous study of the 90° and 180° domain 
wall, we provide strong support, at least in the case of 
defect-free walls at low temperature, that both domain 
walls are of comparable width. However, we find that 
the barrier for the motion of the 90° domain wall is ex- 
tremely small, suggesting that at non-zero temperatures 
the domain walls may fluctuate strongly. 

Many experimental techniques are capable of reveak 
ing the domain microstructure of a ferroelectric crystal,u 
but it is much more difficult to image the domain bound- 
aries and determine their structure with atomic-level ac- 
curacy. Most studies so far have used high-resolution 
transmission electron microscopji-ZHRTEM) to obtain di- 
rect images of the domain wallsJj^llil In a detailed study 
of 90° donw,in walls in BaTiOa and Pb(Zro.52Tio.48)03 
Tsai et aln observed a diffuse dark contrast about 4- 
10 unit cells wide in their HRTEM images of domain 
boundaries, which they attributed to deviations from the 
Bragg condition through lattice distortions and ionic dis- 
placements. But no further analysis of these images was 
done. Combining HRTEM and quantitative image anal- 
ysis, Stemmer et al.B measured the width of the 90° do- 
main wall in PbTiOa to be 10±3 A, and they gave an es- 
timate for the domain- wall energy of 50 mJ/m^. This ev- 
idence for a very narrow 90° domain wall in PbTiOa was 
further confirmed by Foeth et alJid From their HRTEM 
images they obtained a width of 15±3 A at room temper- 
ature, and from the analysis of the thickness fringes of 
weak-beam transmission electron microscopy (WBTEMl 
they found a value of 21±7 A. In contrast, Floquet et al.til 
concluded, based on their combined study of HRTEM 
images and a careful analysis of the X-ray diffraction 
pattern (XRPD) of BaTiOa, that 90° domain walls are 
regions of width 40-60 A where the crystallographic dis- 
continuity is accommodated by irregular atomic displace- 
ments. For micron-sized BaTiOa grains they even pro- 
posed a wall thickness of 140 A. 

This paper is organized as follows. In Section |l| we de- 
scribe the technical details of our computational method 
and the geometry of the supercells used to model the do- 
main walls. In Section III and IV we present our results 
on the 180° and 90° domain walls, respectively. Finally, 
the paper concludes with a summary. Section ^ in which 
we also discuss the results in the context of the previous 
experimental work. 



II. THEORETICAL DETAILS 

A. Method of calculation 

For our calculations we have used a first-principles 
plane-wave pseudopotential method based on density- 
functional theory (DETjEj within the local-density ap- 
proximation (LDA).t3 Ab-initio simulations of this kind 
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a \A] 


3.86 


3.90 


c[A] 


4.04 


4.15 




0.71 


0.72 


C.(Ti) 


0.34 


0.33 


UOi) 


-0.39 


-0.35 


UOii) 


-0.39 


-0.35 


UOiu) 


-0.27 


-0.35 


ICI [A] 


0.38 


0.43 



TABLE I. Computed and experimental values of the struc- 
tural parameters for bulk PbTiOa in the tetragonal phase, a 
and c are the lattice constants, and the atomic displacements 
relative to the cubic positions are given as a unit vector ^ 
times an amplitude |^|. For the labels of the oxygen atoms 
see Fig. Q(a). 
"Ref. m 



have be applied successfully many times to explaia. 
the origin of ferroelectricity in perovskite compounds ,Efl 
to determine ground-state structureal^ and to repro- 
duce sequences of phase transitions.E£lrTke spontaneous 
polarizatiop, LZl Born effective charges,ta't2l |piezoelectric 
coefficient^^! and other material parametersEil have been 
calculated in excellent-affljeement with experiment. Also 
defectsc3 and surfacesEaEJ in ferroelectric materials have 
been studied using DFT methods. 

The pseudopotentials were of the Vanderbilt ultrasoft 
typecj with the semicore Pb 5d and Ti 3s and 3p states 
explicitly treated as valence states (in contrast to Ref. o). 
A conjugate gradient technique as described in Ref. ]l5| 
was employed to minimize the Hohenberg-Kohn total en- 
ergy functional. For the plane- wave expansion of the elec- 
tron wavefunctions a cutoff energy of 25 Ry was used. 
The same pseudopotentials have already been used in 
several previous studies-ou-PbTiOa where their accuracy 
has been demonstrated,^^^^ and where it has been shown 
that a cutoff energy of 25 Ry is sufficient to obtain very 
well converged results. The resulting computed struc- 
tural parameters for the tetragonal phase of PbTiOa are 
given in Table ||. The lattice constants a and c are under- 
estimated by 1 % and 2.6%, respectively, which is typi- 
cal for LDA. The deviations in the atomic displacements 
of the oxygen ions between theory and experiment are 
mainly due to the fact that the X-ray data were ana- 
lyzed under the assumption that the displacement is the 
same for all three oxygen atoms. 

The construction of appropriate supercells for the 
study of the domain walls of interest will be detailed 
in the following subsection. All atomic configurations 
were fully relaxed by minipiizing the atomic forces using 
a variable-metric scheme.EZI Convergence was assumed 
when the forces on the ions were less than 0.005 eV/A. 
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B. Domain wall geometries 

PbTiOa belongs to the important group of ferroelectric 
materials based on the cubic perovskite structure. At the 
temperature of 765 K PbTiOs undergoes a single phase 
transition from the paraelectric cubic phase to a tetrago- 
nally distorted ferroelectric phase, Fig. |^(a), and remains 
in this tetragonal phase down to T= 0. In this state, six 
energetically equivalent orientations of the spontaneous 
polarization exist. The domain walls between these vari- 
ants can be regarded as twin boundaries on low-index 
lattice planes. With the additional constraint that the 
normal component of the polarization should be contin- 
uous across the domain wall, so that no net interface 
charge is present, there are two allowed types: (i) twins 
on (100) planes with parallel polarization of opposite ori- 
entation in the neighboring domains (180° domain wall); 
and (ii) twins on (101) planes with the polarizations on 
either side of the domain wall being almost perpendicular 
to one another (90° domain wall). These situations are 
illustrated in Figs, ^b) and (c), respectively. 

For 180° domain walls, two possible high-symmetry 
cases exist. These are Pb-centered and Ti-centered do- 
main walls, resulting from twinning on PbO- and Ti02- 
planes, respectively. In each case one of the metal cations 
acts as a center of inversion symmetry. Of course, we 
have to check whether one of these symmetries is actu- 
ally present, or whether the atomic relaxation leads to a 
lower-symmetry structure; we will return to this point in 
Sec. m. 

For the 90° domain wall, we can imagine obtaining 
a starting structure by twinning on either Pb-Ti-0 or 
0-0 types of (101) planes. However, in this case the 
"flow" of the polarization establishes an intrinsic differ- 
ence between the "upstream" and "downstream" sides 
of the domain wall (left and right sides, respectively, in 
Fig. 0(c)). Thus, the two sides of the domain wall can- 
not be related by any symmetry operation, and there 
is no possible scenario in which the position of the do- 
main wall is determined by symmetry. In other words, 
we expect that if we start from either starting guess, the 
structure will relax to some ground-state structure lo- 
cated at some intermediate position between the Pb~Ti- 
O and 0-0 limits. There is thus no sharp distinction 
between a Pb-Ti-0 or an 0-0 centered domain wall in 
the 90° case. In our calculations, we start the relaxation 
from the structure obtained by twinning on the Ti-Pb- 
O plane. (As we shall see, the relaxed structure ends up 
being centered much closer to the 0-0 plane, so clearly 
we introduce no bias by doing so.) 

Since periodic boundary conditions are used in plane- 
wave pseudopotential calculations, we cannot study sin- 
gle domain walls in isolation. Instead we have to build 
supercells containing a domain structure that can be re- 
peated periodically in three dimensions. As a conse- 
quence, we always have to include two domain walls in a 
supercell. 



(a) 





Pb • Ti O O 



(b) 



DW 
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(c) 




DW 
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FIG. 1. (a) Unit cell of the cubic perovskite structure and 
its projection along the [010] direction; arrows indicate atomic 
displacements in the ferroelectric tetragonal phase, (b-c) Su- 
percell geometries containing 8 perovskite primitive cells for 
(b) 180° and (c) 90° domain walls (DW's); atoms are omit- 
ted and only solid lines connecting the Pb atoms are drawn. 
Supercell boundaries are indicated by dashed lines. 



For the simulation of the 180° domain wall, we have 
used supercells consisting of iV x 1 x 1 perovskite unit cells 
stacked in the x direction and containing one "up" and 
one "down" domain, each N/2 lattice constants wide, as 
shown in Fig. |l](b) . We take the domain wall to lie in the 
y-z plane and the polarization to point in the z-direction. 
Supercells with iV=6, 8 and 10 were employed, contain- 
ing 30, 40 and 50 atoms, respectively. The dimensions of 
the supercell were kept fixed at the theoretical equilib- 
rium lattice constants computed for the bulk tetragonal 
phase (see Table ||). For all Brillouin-zcme integrations 
a (2,4,4) Monkhorst-Pack k-point meshca was used; this 
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k-point mesi-has previously been shown to be sufficient 
for PbTiOa.y 

The geometry of the supercell containing two 90° do- 
main walls is sketched in Fig. |l|(c). The tetragonal shape 
of the PbTiOa unit cell causes the angle between the c 
axes on either side of the domain wall to deviate slightly 
from 90°. Instead, to form a coherent twin interface, 
an angle of 2arctan(a/c) is geometrically required. This 
gives rise to a very characteristic bending of the [100] 
and [001] atomic rows (a- and c-rows) across the do- 
main wall; in fact, this bending is precisely the signature 
through, which 90° domain walls are identified in HRTEM 
imagesJj Again, we have taken the domain wall to lie in 
the y-z plane and the x direction to be perpendicular to 
the wall. The supercell is orthorhombic with dimensions 
of Nc/ yjl + {c/aY, a, and a-y/1 + {c/aY in the x-, y-, 
and z-directions, respectively. We used supercells con- 
sisting of N—Q, 8, 10, 12 and 14 perovskite unit cells, 
thus containing up to 70 atoms. The theoretical LDA 
values from Table || were taken for the lattice parameters 
a and c. For one supercell we also relaxed the lattice con- 
stant in the x direction, but the resulting expansion of the 
supercell in the x direction was very small (less than 0.8 % 
of a), and the changes in the atomic displacements and in 
the domain- wall energy were negligible. Two symmetry 
operations were imposed during the atomic relaxation: a 
mirror operation across the y=0 plane, and a glide mirror 
operation across the z—Q plane with a translation of half 
the superlattice dimension in the x direction. The first 
symmetry corresponds to a physical assumption that the 
polarization does not rotate out of the y=Q plane, while 
the second is just a technical convenience to insure that 
both domains in the supercell stay equivalent. For the 
Brillouin zone-sampling again a (2,4,4) Monkhorst-Pack 
k-point meshcj was used. 



III. RESULTS FOR THE 180° DOMAIN WALL 

We have investigated both the Pb- and Ti-centered 
domain walls. The relaxation of the atomic structure 
was started from the ideal twinned configurations, and 
the inversion symmetry centered on a cation in the do- 
main wall was enforced during the structural optimiza- 
tion process. This prevents a Ti-centered domain wall 
from transforming into a Pb-centered wall or vice versa. 



A. Domain wall energy 

The results for the domain-wall energies from our vari- 
ous supercells are summarized in Table ||. The reference 
energy of the bulk structure (without domain walls) was 
calculated from a single perovskite unit cell using a k- 
point mesh equivalent to the one used for the supercell. 
Values for BaTiOs are given for comparison; these were 
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132 


169 


7.2 


16.5 
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132 


169 


7.4 


16.7 


10 


132 


169 


7.5 


16.8 



TABLE II. Calculated 180° domain-wall energies 
(mJ/m^) using supercells of A'' perovskite primitive cells. 



calculated by the same procedure except that a (2,6,6) 
k-point mesh was employed. 

All relaxed structures have inversion symmetry about 
a cation in the domain wall (as well as mirror symmetry 
about the y=Q plane). The domain- wall energy is low- 
est for the Pb-centered domain wall; this is therefore the 
preferred wall configuration. The same ordering holds for 
BaTiOa, but the domain- wall energies are significantly 
lower. As can be seen from the Table, the domain-wall 
energy is well converged even with the smallest separa- 
tions of the domain walls. |-| 

In contrast, Poykko and ChadiQ found that domains 
have to be at least four lattice constants wide before 
the up/down domain pattern becomes stable. Further- 
more, their calculations indicated that the Pb-centered 
domain-wall energy could be lowered by almost a factor 
of two by breaking the inversion symmetry about a Pb 
atom in the domain wall. In their lowest-energy struc- 
ture, in which there are large relaxations of Pb and O po- 
sitions in the domain- wall plane, the only remaining sym- 
metry was a 180° rotation about a z-axis (C2z) passing 
through the center of a domain. In view of those results, 
we made careful tests in which we lowered the symmetry 
of our supercell in two ways: (i) as in the final result of 
Ref. ^, we kept only the C2z symmetry about the center 
of a domain; and (ii) as in the case of the 90° domain 
wall, we kept only the glide mirror plane at z—Q with a 
half-supercell translation along x. We then distorted the 
atomic configuration in several ways compatible with the 
lowered symmetry and started a new atomic relaxation. 
In all cases, the structure relaxed back toward the high- 
symmetry centrosymmetric structure. The relaxation is 
very slow, indicating a very flat energy surface, but we 
never found a lower-energy structure with reduced sym- 
metry. 

Regarding the domain-wall energy, our value of 
132mJ/m^ for the Pb-centered wall is much lower than 
the value of 270mJ/m^ obtained in Ref. ||when inversion 
symmetry is imposed. Strangely, it is much closer to the 
150mJ/m^ value reported in Ref. ||when inversion sym- 
metry breaking is allowed. For the Ti-centered domain 
wall our wall energy of 169mJ/m^ again lies below their 
value of 220mJ/m^. The case of a Ti-centered wall with 
broken inversion symmetry cannot be considered, since 
relaxation of the domain wall would then just tra nsform 
it into a Pb-centered wall, as we will show in Sec. Ill C 
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FIG. 2. Change of polarization across the (a) 

Pb-centered, and (b) Ti-centered, 180° domain wall. Rz de- 
notes the ferroelectric distortion of each lattice plane in the z 
direction, in units of the distortion associated with the bulk 
spontaneous polarization. 
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B. Atomistic domain wall structure 

We begin by analyzing how much the ferroelectric dis- 
tortions along the z direction are changed by the presence 
of a domain wall. To give a quantitative description we 
compute an average ferroelectric distortion i5fe for each 
atomic plane parallel to the domain wall. As a measure 
for (5fe we take the displacement of the metal atonujeU. 
ative to an oxygen atom in the PbO or Ti02 plane.Ej'0 
(In the case of a Ti02 plane, we choose the oxygen ion 
lying along the z direction from the Ti atom.) We then 
define the parameter Rz as the ratio between (5fe of the 
specified lattice plane and its value in the undistorted 
ferroelectric bulk phase. Rz is a very sensitive indicator 
of the ferroelectric order: Rz is zero as long as the atoms 
are in their cubic positions, and tends to unity as the full 
bulk ferroelectric distortion is attained. 

The results for the distortion Rz are illustrated in 
Fig. H for both Pb-centered and Ti-centered domain 
walls. Only the results from the 50-atom supercells are 
shown; results for smaller cells are essentially identical. 
We focus our discussion here on the more physical Pb- 
centered domain wall. We find that Rz— 0.80 already for 
the Ti02 first-neighbor plane (to be compared with the 
value of 0.73 reported in Ref. ^), and the ferroelectric dis- 
tortion is essentially fully recovered to its bulk value by 
the PbO second-layer plane. The orientation of the polar- 
ization thus changes abruptly over a distance of less than 
two lattice constants, leading to a very narrow domain 
wall. The narrowness of the 180° domain wall had earlier 



FIG. 3. The z-displacement of the atoms relative to the 
Pb atom in a Pb-centered 180° domain wall as function of 
the distance to the wall, calculated with different supercell 
sizes A''. Integer distances represent PbO planes, half-integer 
distances Ti02 planes. Labels of oxygen atoms correspond to 
those in Fig. ^(a). 

been predicted based on phenomenological modelsEJ and 
confirmed using atomic force microscopy (AFM).c3 The 
results for the Ti-centered wall are qualitatively similar. 

In the next step we analyze for the Pb-centered domain 
wall how the atomic rows are aligned across the wall. As 
already indicated in Fig. |^, there is a considerable offset 
between [100] atomic rows to the left and right of the 
domain wall. Fig. ^ shows the detailed results for the z 
displacements of the atom rows relative to the Pb ion in 
the domain wall. The discontinuity in the [100] rows is 
largest for the Pb atoms, where we find a jump of 14.5 % 
of the lattice constant c. This offset of 0.6 A should be 
clearly visible in an HRTEM image, thus allowing identi- 
fication of 180° domain walls having polarization parallel 
to the surface in the neighboring domains. 

It is not surprising that the largest offset occurs for the 
Pb rows, since the Pb displacement contributes most to 
the ferroelectric soft- mode vector (see Table ||) . By com- 
pensating for the different contributions of the atoms to 
the ferroelectric soft-mode vector, we can define a kind of 
"geometrical offset" between neighboring domains. We 
do this by going to the center of each domain and sub- 
tracting the bulk ferroelectric displacement pattern from 



Q 



170 




1 1 I I 1 



1 ^ — ^ 


160 




X / ■ 




(132mJ/m') 


j/ Ti-centered 


150 




domain wall 


- Pb-centered 






domain wall ^ 
1 / 


(169mJ/m') 


140 


- fx 




130 







iV 


6 


8 


10 


12 


14 


-Bdw 


29.4 


32.5 


34.6 


34.7 


35.2 


-f^barricr 


0.6 


1.1 


1.4 


1.5 


1.6 



0.5 1 

domain wall position X 

FIG. 4. Energy profile along the estimated path for the 
motion of a 180° domain waU, calculated with supercells of 
10 perovskite unit cells. 



the atomic coordinates, thereby shifting all atoms into 
essentially cubic positions. The offsets of the resulting 
atomic rows are now the same for all five atoms, and 
this joint offset therefore provides a reasonable measure 
of the geometrical misalignment. It turns out that this 
geometrical offset is very small; we found values of 0.010 c 
and 0.008 c for Pb-centered and Ti-centered walls respec- 
tively. 



C. Barrier for domain wall motion 

Many models have been proposed in the literature for 
the motion of domain walls in crystals. Most realistic 
models involve the formation and propagation of kinks of 
various shapes, allowing the domain wall to shift gradu- 
ally to its new position. It is not our aim here to comment 
on or distinguish between these models; instead, we con- 
centrate only on extracting one parameter that may be 
needed in many models of this type. Specifically, we com- 
pute the barrier height for coherent motion of the entire 
domain wall to its neighboring lattice position. Thus, our 
calculations always preserve full translational symmetry 
parallel to the domain wall. 

We made a simple estimate for the path in configura- 
tion space along which a Pb-centered domain wall may 
gradually transform into a Ti-centered wall as follows. 
We took the atomic configurations ^Pb and ^xi of Pb- 
and Ti-centered walls on neighboring lattice planes and 
formed the linear interpolation ^pb + A(^Ti— ^Pb), where 
A has the interpretation of being a reaction coordinate 
along the path. The energy profile along this path is 
shown in Fig. ^. The configurations of the Pb- and Ti- 
centered domain wall are local minima and maxima of 
the energy curve, respectively, with no barrier between 
them. The Ti-centered wall is therefore metastable and 
will transform spontaneously into a Pb-centered domain 



TABLE III. Calculated 90° domain-wall energies _Edw 
and barrier height for domain wall motion -Ebarrier (mj/m^) 
as function of supercell size A'^. 



wall. In our calculation the Ti-centered wall could only 
be stabilized by imposing a center of inversion symmetry 
in the domain wall. Moreover, Fig. ^does not give any in- 
dication that the energy of the Pb-centered domain wall 
can be lowered by breaking the inversion symmetry. Hav- 
ing identified the Ti-centered domain wall as the saddle 
point configuration for the motion of a 180° domain wall, 
the barrier height for a jump to the nearest-neighbor lat- 
tice position is given by the difference between the wall 
energies of the Ti- and Pb-centered walls. We find this 
difference to be 37mJ/m^. 



IV. RESULTS FOR THE 90° DOMAIN WALL 

A. Domain wall energy 

For an accurate calculation of the domain-wall energy, 
it is indispensable to reduce all effects of systematic errors 
as far as possible, since an uncertainty of only lOmeV in 
the energy difference between the supercell calculation 
and the reference structure already leads to an inaccu- 
racy of 7.5mJ/m^ in the domain- wall energy. Because of 
the particular shape of the superce lls re presenting the 90° 
domain wall (as described in Sec. [IB), it is not possible 
to construct equivalent k-point meshes for the supercells 
and single perovskite unit cells (in contrast to the 180° 
case) . To get good energy differences between the relaxed 
structures and the bulk reference, we applied the follow- 
ing procedure for calculating the domain- wall energies. 
First, for the full atomic relaxation of the supercells, we 
used our regular (2,4,4) Monkhorst-Pack k-point mesh. 
Then, to reduce the k-point errors further, we reran the 
total-energy calculation of the supercell with a denser 
(2,6,6) Monkhorst-Pack mesh. The reference energy of 
the bulk structure was computed from a single perovskite 
unit cell using a k-point mesh with the same density and 
an orientation that corresponds as closely as possible to 
the supercell (2,6,6) Monkhorst-Pack mesh. Specific.ally, 
we employed a (TV, 6, 6) Moreno- Soler k-point meshEJ ori- 
ented along the [101], [010], and [101] directions (see 
Fig. 0(c)). 

The results for the 90° domain-wall energies from the 
various supercells are given in Table III, The conver- 
gence of the wall energy with supercell size is somewhat 
slower than for the 180° case, but the energy is well 
converged for a wall separation of five perovskite unit 
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cells. The calculated value of 35mJ/m^ is in reason- 
able agreement with a wugh estimate of 50 mJ / m^ from 
a HRTEM experiment.tl The energy of the 90° domain 
wall is thus calculated to be about four times lower than 
that of its 180° counterpart, and we can conclude that 
the 90° domain wall is the most stable wall configuration 
in PbTiOs. 

B. Domain wall structure and polarization profile 

In the case of the 180° domain wall the analysis of 
how the ferroelectric distortions change in the vicinity of 
the domain wall was relatively straightforward, since the 
question could be treated as a one-dimensional problem 
(the relaxation of the atoms in the x-direction, perpen- 
dicular to the domain wall, being completely negligible). 
In contrast, for the 90° domain wall the two-dimensional 
relaxations of the atoms lead to a complex relaxation pat- 
tern for which an analysis in terms of average ferroelectric 
distortions (5fe and their ratio relative to the bulk value 
Rz would be difhcult to repeat. 

Therefore, to simplify the discussion of the changes in 
the ferroelectric distortions, we have chosen to visualize 
the atomic displacements by translating them into a "po- 
larization P*^*-' per unit cell i" as follows. First, we have 
to decide what we mean by a "unit cell" in the context of 
the domain- wall structure. Figure || shows three possible 
ways of doing this. For the 'Ti-centered' choice, for exam- 
ple, the unit cell contains the central Ti atom, but the six 
neighboring O atoms are each shared by two unit cells, 
and the eight neighboring Pb atoms are each shared by 
eight unit cells. Thus we assign weights WTi=l, ?i'o=l/2, 
and wph—l/8 to describe the association of atoms to 
the unit cell. Two other choices are the 'Pb-centered' 
and 'M-centered' (metal-centered) cells, also illustrated 
in Fig. ||. The M-centered cell is the most compact in the 
X direction; it therefore gives the finest possible resolu- 
tion to the function P(a;) in the case of the 90° domain 
wall, which consists of alternate Pb-Ti-0- and 0-0- 
type (101) planes stacked along x (see Fig. |l|(c)). Next, 
we attach a polarization to each unit cell i via the assign- 
ment 

p^'' = E«^"Z:.uW . (1) 

Here e is the electron charge, fie the volume of the five- 
atom bulk unit cell, index a runs over all atoms of unit 
cell i, Uq is the displacement of atom a from its cubic 
position, and Wa are the weight factors introduced above. 
The Z* are the Born effective charge tensors, which are 
approximated by their cubic bulk values as discussed be- 
low. 

For the case of a homogeneous system having the prim- 
itive five-atom periodicity, Eq. (Q) is independent of cell 
index i, and even of the choice of unit cell in Fig. g; in 



Ti-centered Pb-centered M-centered 




FIG. 5. Sketch of the three different choices of unit cells 
used for the calculation of local polarization via Born effec- 
tive charges. Shown is the projection on the a;-z-plane (see 
Fig. 0(a)). 

this case, Eq. (|l]) simply gives the bulk macroscopic po- 
larization. In the inhomogeneous case the spatial decom- 
position indicated by Eq. (|l|) is somewhat arbitrary, but 
we think that it is a rather natural choice. Of course, 
the cell-by-cell polarization now depends on the choice 
of unit cell in Fig. ^ Moreover, in principal it would be 
necessary to calculate the Born effective charge tensor for 
each atom in the neighborhood of the domain wall, since 
these will differ from the bulk-like values. Here, we make 
the further assumption of using the cubic bulk effective 
charge tensors in Eq. (|l|) . Previous work has shown that 
the Z* matrix elements change by less than 10 % inpftoing 
from the cubic to the tetragonal phase of PbTiOa.ta We 
expect a variation of similar magnitude for the atoms in 
the interior of the domain wall. Moreover, the sum rule 
for the Born effective charges guarantees that a simple 
uniform translation of the entire contents of the domain- 
wall supercell will not change the P'*-' values at all. We 
therefore think this is a reasonable approximation. How- 
ever, it is clear that the results presented on the basis of 
this approximation should be regarded mainly as a visu- 
alization of the complex atomic relaxation pattern, not 
as a realistic calculation of the spatial polarization field. 
In fact, one should keep in mind that defining a local po- 
larization is an intrinsically subtle task. The polarization 
is a macroscopic quantity, and there is always some am- 
biguity when this concept is transferred to an atomistic 
level. 

For the Born effective charges Z* we took the theoreti- 
cal values for the cubic structure of PbTiDa as calculated 
by Zhong, King-Smith and Vanderbilt.ll3 Using Eq. (|l|) 
together with the atomic displacements from Table | we 
get a bulk polarization of 89.5/iC/cm^. This result is 
only slightly lacaer than the value of 81.2 /xC/cni2^ cal- 
culated directlylHil with the Berry-phase approach.!!!! All 
three unit cells in Fig. || are centered around a Pb-Ti-0 
plane; we therefore assign the calculated local polariza- 
tions to the corresponding Pb-Ti-0 planes and plot the 
polarization profiles as a function of the positions of these 
planes. 

Such a plot is shown in Fig. ^ The corresponding 
curves are least-square fits to functions of the form 
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;dw 



and 



Pz{x) = Pj^o tanh 



X — Xq 



;dw 



(2) 



(3) 



where Px,o and Pz,o specify the polarization deep in 
the domain interior, and xq and 2^dw correspond to 
the center-position and width of the domain wall re- 
spectively. (The factor of two in 2^dw is a convention 
that is included to facilitate comparisop-with experimen- 
tal reports that use this convention. BO) Eq. (^) is the 
soliton solution for a domain wall in a one-dimensional 
fourth-order Landau-Ginzburg-type theory.Q However, 
since PbTiOs is a proper ferroelectric (or improper fer- 
roelastic), it would be more appropriate tO|-jiS€ the result 
of a sixth-order Landau-Ginzburg theory.uc2l The solu- 
tion for a domain wall would then be 



Pzix) = P, 



sinh 



/ X-Xq \ 
\ C DW / 



B + sinh- 



I X-Xq \ 
\ C DW / 



1/2 



(4) 



which is identical to Eq. (p),for B—\. Since B—IA has 
been reported for PbTiO^JEI the deviation between the 
two expressions is not very large. 

In Fig. ^ the component of the polarization parallel to 
the domain wall P^ changes its sign very abruptly over a 
distance of less than three lattice spacings between Pb- 
Ti-O-type (101) planes. Furthermore, it can be seen that 

(i) 

the results for the local polarizations Pz do not depend 
strongly on the choice of unit cell. For the Ti- and Pb~ 
centered cell they are almost indistinguishable, and for 
the M~centered cell the polarization changes even more 
rapidly. Thus, as expected, the M-centered cell gives 
the finest resolution for the polarization rotation in the 
domain wall. 

From the fits to Eq. (^ we find in all three cases 
a domain wall position that is very nearly halfway be- 
tween two Pb-Ti-0 planes. Thus, from an electrostatic 
point of view, the domain wall can be regarded as be- 
ing essentially centered on an 0-0-plane. For the do- 
main wall width we find a value of 2^ dw = 5±0.5 A. This 
even somewhat narrower than the results, p£ 10±3 A and 
15±3A suggested by recent experiment .B'tll However, it 
has to be kept in mind that the experiments were done at 
room temperature, whereas our calculations are for zero 
temperature. 

Figure ^ shows that the component of the polariza- 
tion perpendicular to the domain wall is almost constant, 
showing only a small reduction as it crosses the domain 
wall. The local polarizations Px scatter a little bit more 
with the choice of unit cell than in the case of the par- 
allel components Pz^\ and it is more difficult to obtain 
meaningful fits to analytic curves. (Note also that the 
changes of Px occur on a much finer scale than those 
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FIG. 6. Polarization profile across a 90° domain wall, cal- 
culated with a supercell of 14 perovskite unit cells. The sym- 
bols represent the local polarizations Pi*' and Pi'' as obtained 
using the three different choices of the unit cell illustrated in 
Fig. ^ solid, dotted and dashed lines are fits using Ti-, Pb- 
and M-centered unit cells, respectively. 

of Pz^\) Despite the larger variations between the fit- 
ted curves, in all three cases the minimum of Pi*^ lies 
near the same 0-0 plane where Pz crosses through zero. 
Thus, the analysis of the Px curves reinforces the in- 
terpretation that the domain wall is centered on an 0-0 
atomic plane. 

In a full three-dimensional treiatment of a Landau- 
Ginzburg model Cao and CrossB found a quasi-one- 
dimensional solution for the polarization profile across 
90° domain walls with a constant perpendicular compo- 
nent Px{x). However, this analytic solution was obtained 
for a special choice of parameters in the Landau-Ginzburg 
functional, and our result of a varying Px{x) is not in con- 
tradiction with their final system of differential equations 
for the general case. 

Finally, we note that while the reduction in Pi'^ in the 
interior of the domain wall appears rather small, it may 
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FIG. 7. Energy profile along the estimated path for the 
motion of a 90° domain wall, calculated with supercells of 10 
perovskite unit cells. 



have important consequences for the electrost atics o f 90° 
domain walls. We return to this point in Sec. IV D. 



C. Barrier for domain wall motion 

As in the case of the 180° domain wall, we address 
here only the situation in which the entire 90° domain 
wall moves coherently from one (101) lattice plane to the 
next, preserving periodicity parallel to the wall. For the 
path in configuration space along which this transfor- 
mation occurs, we again formed the linear interpolation 
in + — £,n) between the atomic configurations <^„ 

and ^n+i of two fully relaxed supercells with domain wall 
positions on nearest-neighbor (101) lattice planes. 

The energies of these atomic configurations are plotted 
as a function of A in Fig. 0. The energy profile is found to 
be a rather featureless cosine-like curve with a maximum 
at A«0.5. The energy at this maximum provides an es- 
timate of the barrier height iJbarricr for the translation 
of a 90° domain wall to a neighboring lattice plane. To 
be more precise, it gives an upper bound for the actual 
barrier, since our assumed path is only an approxima- 
tion to the true one. To refine this estimate further, we 
inspected the forces remaining on the atoms at our esti- 
mated saddle-point configuration at A = 0.5. These forces 
were very small, only slightly above the convergence cri- 
teria used in the structural optimization process. Fol- 
lowing a line minimization in configuration space along 
the direction given by this force vector, the barrier height 
was lowered by less than 0.15 mJ/m^. 

The results for i^barricr obtained using various super- 

At about 1.6mJ/m^ 



cells are listed in Table ID 



this 

barrier is extremely small. It is only ~4 % of the cre- 
ation energy of the 90° domain wall, to be compared with 
the 180° case where the barrier height is ~30 % of the 
domain-wall energy. Or, to put it in relation to thermal 



fluctuations, a barrier height of 1.6mJ/m^ over the area 
of 12 perovskite unit cells corresponds roughly to fc^T 
at room temperature. Therefore, as long as no pinning 
centers are present, the domain wall may be expected 
to fluctuate rather freely except at very low tempera- 
ture. Clearly, it will be important to take this behavior 
into account when interpreting experimental probes of 
the widths and other properties of the 90° domain walls. 



D. 



Electrostatic potential step across the domain 
wall 



In Sec. II B we pointed out that the left and right sides 
of a 90° domain wall are not related by any symmetry 
operation. It follows that the charge distribution across 
the domain wall may have the character of a dipole layer, 
such that the macroscopic electrostatic potential would 
experience a step across the domain wall. As a general 
rule, if something is not prevented by symmetry, it should 
be present. For example, all wurzite crystals have some 
spontaneous polarization, since it is not prevented by 
symmetry; the only question is how large it may be. By 
the same token, since a step in the electrostatic potential 
is not ruled out by symmetry, it should exist for any 90° 
domain wall in a ferroelectric material. The issue then 
is to determine the sign and magnitude of such a step. 
Of course, it might have turned out to be insignificantly 
small, but we s how now that this is not the case. 

In Sec. IV B| we found that the component of the 
polarization perpendicular to the domain wall is slightly 
reduced in the interior of the domain wall. A decreasing 
(increasing) polarization on the left (right) side of the do- 
main wall leads to an accumulation of positive (negative) 
charge, as is sketched in Fig. ||(a). Thus, the dip in 
visible in Fig. |^(a) is a direct fingerprint of the existence 
of such a dipole layer. As a consequence, the electrostatic 
potential should jump from one side of the domain wall 
to the other, leading to band offsets for both the valence 
and conduction bands (see Fig. ^(b)). Mathematically, if 
v{r) is the electron potential energy, Poisson's equation 



V^i;(r) = -47re divP(r) 
can be manipulated to the form 

(Pv/dx^ = -4:TTedP^/da 



(5) 



(6) 



where barred quantities are y-z planar averages, e.g., 

v{x) ^ ^ J^v{r) dydz . (7) 

Then it follows that the offset Av in the electrostatic 
potential for the electrons is 



Av = 



-47re J {Px{x)-Px,o) dx 



(8) 
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FIG. 8. Schematic illustration of a periodic array of 90° 
domain walls (DW's). (a) The component of the polar- 
ization perpendicular to the DW, and the resulting induced 
charge and electron potential v(x), assuming zero macroscopic 
electric field in the interior of each domain, (b) Corresponding 
band offsets, (c) Potential vix) in the case of the correspond- 
ing supercell calculation using periodic boundary conditions, 
i.e., in which the supercell-averaged electric field is zero. 

Using the fits to Eq. ^ illustrated in Fig. ^(a), we find 
Az; = 0.18eV. 

The dipole barrier can be estimated in a completely 
independent way by directly inspecting the electrostatic 
potential v{v) that is calculated in the Kohn-Sham su- 
percell calculation. We first compute planar averages as 
in Eq. and then take sliding-window averages 



v[x) 



L 



L/2 



v{x') dx' 



(9) 



-L/2 



over a window whose width is the distance L between 
two (101) planes. The resulting v{x) is shown in Fig. ^. 

In order to interpret such a figure, it is important 
to understand the role of periodic boundary conditions 
in the supercell calculation. The situation sketched in 
Fig. ^(a) is incompatible with such boundary conditions; 
although the charge density is periodic, the potential is 
not. In all practical DFT/LDA calculations, the poten- 
tial is obtained by solving Poisson's equation under peri- 
odic boundary conditions, resulting in a periodic poten- 
tial as sketched in Fig. ||(c). In fact, the potentials of 
Fig. ||(a) and (c) differ by just a linear potential; in other 
words, in the latter case an artificial uniform electric field 
has been added in order to cancel the potential jump at 
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FIG. 9. Sliding-window averaged electron potential v{x) 
across the 90° domain wall, calculated with a supercell of 12 
perovskite unit cells. 



the supercell boundary. Thus, instead of a step-like be- 
havior as in Figs. ||(a), we expect a zig-zag behavior as 
shown in Fig. [^(c). 

Despite using the sliding-window average, the electro- 
static potential in Fig. || still contains substantial noise 
arising from the different atomic relaxations in the atomic 
planes parallel to the domain wall. Nevertheless, the ex- 
pected zig-zag behavior is evident, as indicated by the 
dashed lines. In particular, the potential jump is very 
close to the expected location where Pz changes sign be- 
tween two Pb-Ti-0 planes (compare with Figs. ^ and||), 
and the height Aw = 0.16 eV of the jump is in very good 
agreement with our previous estimate of 0.18 eV based 
solely upon atomic positions, without any information 
about the electrostatic potential. 

We end this section with a few remarks on the effects 
of the artificial electric field that is introduced in our 
supercell calculation by the dipole layer at the domain 
wall. The strength of this field depends on the supercell 
size and decreases only linearly with reciprocal domain- 
wall separation. It will therefore have some effects on the 
atomic relaxations. For example, if we look at the local 

we 



polarization in the center of the domains. Table IV 



see that the angle a of the polarization with the x axis is 
not arctan(c/a) — 46.3°, as would be expected for an ide- 
ally twinned crystal, but is slightly reduced. The tilt of 
the polarization toward the x axis can be seen to decrease 
with supercell size, and the angle a converges toward it 
ideal value. Furthermore, the amplitude of the polariza- 
tion in the center of the domains is slightly larger than 
in the undisturbed bulk system. It is natural to inter- 
pret this excess polarization as having been induced by 
the artificial electric field. A rough estimate shows that 
for the electric field that is present, such an increase in 
the polarization can be expected if we assume a dielec- 
tric constant of roughly 50, which is not unreasonable. 
Overall, these small effects of the artificial electric field 
are most likely the origin of the slower convergence with 
supercell size of the domain-wall energy and migration 
barrier height compared to the 180° case, where no such 
electric field is present. 
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6 


74.7 


42.9 


86.1 


29.8 


8 


72.4 


52.4 


89.3 


35.9 


10 


70.8 


57.8 


91.5 


39.2 


12 


69.7 


59.4 


91.6 


40.4 


14 


68.9 


60.4 


91.6 


41.3 


bulk 


61.9 


64.7 


89.5 


46.3 



TABLE IV. The x and z components of the polarization 
(/iC/cm^) in the centers of the domains for different supercell 
sizes TV. P is the norm of the polarization vector and a its 
angle with respect to the x-axis. 



V. SUMMARY 

In summary, wc have studied the properties of 180° 
and 90° domain walls in PbTiOs using density-functional 
theory and the local-density approximation. Metastable 
supercells representing both kinds of domain walls were 
constructed, and the fully relaxed atomic configurations 
of the domain walls were obtained by minimizing the 
residual forces on the ions. 

In the 180° case, the domain wall located on a Pb-0 
lattice plane is lowest in energy, whereas the Ti centered 
wall is a saddle point configuration that spontaneously 
transforms into a Pb-centered wall. The domain wall is 
found to be very narrow, with a width on the order of 
the lattice constant a, and the energy to create a 180° 
domain wall is calculated to be 132mJ/m^, in very good 
agreement with a previous ab-initio study. 

For the 90° domain wall we found a much lower 
domain-wall energy of 35 mJ /m^ . From the polarization 
profile across the domain wall we obtained the position 
of the wall to be halfway between Pb-Ti-0 planes. We 
also find a domain- wall width of 5±0.5 A, the same order 
of magnitude as for the 180° wall. This conclusion about 
the narrowness of the 90° domain wall is in line with two 
recent HRTEM experiments. Together, these works indi- 
cate that the conventional wisdom, which holds that the 
90° domain wall is much broader than its 180° counter- 
part, does not hold true in PbTiOs. 

On the other hand, the barrier height for the coher- 
ent motion of a 90° domain wall to a neighboring lat- 
tice plane is estimated to be very small. The domain 
wall can easily move so that a scenario where the do- 
main wall is pinned at defects and fluctuates widely in 
between is very likely to occur in real samples. The ap- 
parent width of the domain boundary may therefore be 
much broader, raising the question as to what exactly 
we mean by the "width of a domain wall." In this re- 
gard, our result 2^dw = 5±0.5 A for an ideal, defect- free 
wall should be regarded as a intrinsic width of the do- 
main boundary at zero temperature, whereas at finite 
temperatures the domain wall may be smeared out by 
fluctuations. Regarding the sharp HRTEM images of the 



90° domain walls, we also have to take into account that 
the samples in these kinds of experiments are very thin 
(typically 10 nm) so that surface pinning of the domain 
walls could play an important role. 

Finally, we have found an electrostatic dipole moment 
across the 90° domain wall, which causes a step in the 
electrostatic potential and an offset of the valence and 
conduction bands of 0.15-0.2 eV. To our knowledge, the 
existence of such a dipole moment across a ferroelectric 
domain wall has never been reported in literature before. 
It remains a challenge whether such an effect can be ob- 
served experimentally. 
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